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r\l ■ For the anisotropic XY model in transverse magnetic field, we analyze the ground state and its 

concurrence-free point for generic anisotropy, and the time evolution of initial Bell states created 
in a fully polarized background and on the ground state. We find that the pairwise entanglement 

gj , propagates with a velocity proportional to the reduced interaction for all the four Bell states. A 

transmutation from singlet-like to triplet-like states is observed during the propagation. Charac- 
teristic for the anisotropic models is the instantaneous creation of pairwise entanglement from a 

C/3 ^ fully polarized state; furthermore, the propagation of pairwise entanglement is suppressed in favor 

of a creation of different types of entanglement. The "entanglement wave" evolving from a Bell 
state on the ground state turns out to be very localized in space-time. Our findings agree with a 
recently formulated conjecture on entanglement sharing; some results are interpreted in terms of 
this conjecture. 
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I. INTRODUCTION 



A quantum mechanical system possesses additional correlations that do not have a classical 

^ , counterpart. This phenomenon, called entanglement, is probably one of the most astonishing 

^f^ • features of quantum mechanics. Understanding the nature of these non-local correlations has been 

lO [ a central issue in the discussion of the foundation of quantum mechanics. More recently, with 

*/^ . the burst of interest in quantum information processing, entanglement has been identified as an 

'bl ' ingredient for the speed-up in quantum computation and quantum communication protocols^ as 

compared with their classical counterparts. Therefore it is of crucial importance to be able to 

generate, manipulate and detect entangled states. Experimental efforts in this direction have been 

put forward using photons^^, cavity QED systems^, ion traps-, and coupled quantum dots^. Also 

"j^ ' encouraged by the advances in the field of nanoscience, there has been a number of proposals 

Cj I to detect signatures of entanglement. Systems under current study are multiterminal mesoscopic 

I i deviceai^ and Josephson junctions^. 

' O ' Though in many-body systems correlated states naturally appear, a research activity investigating 

C , entanglement in condensed matter systems emerged only recently. Despite conceptual difficulties, 

e.g. in distilling the quantum part of the correlations while the interaction between the subsystems 

is on, many-body systems might become useful for the development of new computation schemes 

and/or communication protocols. As an example we mention the recent proposal by Bose^" to use 

the spin dynamics in Heisenberg rings to transfer quantum states. It is conceiveable that along the 

same lines other quantum information tasks can be implemented as well. 

C^ I An important motivation for us to study the interconnection between condensed matter and 

quantum information is to investigate whether it is possible to gain additional insight in condensed 

matter states from quantum information theoryii. The peculiar aspects of non-local correlations 

become particularly evident when many bodies behave collectively; a prominent example is a system 

close to a quantum phase transitions^, where it was found that entanglement can be classified in 

the framework of scaling theorjiiii^ii^*i&, but also profound differences between non-local quantum 

and classical correlations have been highlighted. In a very recent paper the problem of decoherence 

in a near-critical one-dimensional system was addressedii. The study of entanglement has not 

been devoted only to spin systems, but also to the BCS modeliSiiS, quantum Hall^'^'-l and Boson 

systemsSS. 

In this paper, we examine the evolution of a local excitation bearing entanglement. There are a 
number of questions that can be addressed in this way. In particular we would like to see if there 



is a well defined velocity for the transport of entanglement and how it is related to the propagation 
of the elementary excitations of the spin system. Furthermore: is there a parameter regime which 
favors entanglement transport even over larger distances; what are the time scales for the damping 
of entanglement created initially. Another important question to analyze is how the transport 
of entanglement is influenced by the interference with other entangled states. Finally, we try to 
discriminate pairwise from other types of entanglement. This problem was quantified by Coffman, 
Kundu, and Wootters (CKW) in terms of a conjecture on a measure for residual entanglement^"^. 

The paper is organized as follows: In the next section we will introduce the model Hamiltonian, 
its spectrum and correlation functions needed in the subsequent sections. In section UTTl we present 
the applied entanglement measures. Section Hvl is revisiting the ground state entanglement and the 
entanglement dynamics for the isotropic and anisotropic model is presented in section The final 
section is devoted to conclusions. 



II. THE MODEL 

The system under consideration is a spin- 1/2 ferromagnetic chain with an exchange coupling A 
in a transverse magnetic field of strength h. The Hamiltonian is H = hHs with the dimensionless 
Hamilton operator Ha being 

N N 

if, = -A 5](i + 7)5rsr+i + (1 - 7)5f 5f+i - E ^^r (1) 

1=1 i=l 

where S°' are the spin-1/2 matrices (a — x, y, z) and N is the number of sites. We assume periodic 
boundary conditions. The anisotropy parameter 7 connects the quantum Ising model for 7=1 
with the isotropic XY model for 7 = 0. In the interval < 7 < 1 the model belongs to the Ising 
universality class and for iV = cxd it undergoes a quantum phase transition at the critical value 
Ac = 1. The order parameter is the magnetization in x-direction, (5^), which is different from zero 
for A > 1 and vanishes at and below the transition. On the contrary the magnetization along the 
^-direction, {S^), is different from zero for any value of A. 

This class of models can be diagonalized by means of the Jordan- Wigner transformation^iiSSiSSiSi 
that maps spins to one dimensional spinless fermions with creation and annihilation operators c] and 
C;. It is convenient to use the operators Ai := cj +Ci, Bi := cj — Ci, which fulfill the anticommutation 
rules {Ai, Am} = — {_B;, S,„} — 26im, {Ai,B,n} = 0. In terms of these operators the Jordan- Wigner 
transformation reads S[ = ^AiH^Z^^^ AsB^, Sf = -|B/ni=\^^s, and S[ = -^AiBi. The 
Hamiltonian defined in Eq. |^ is bilinear in the fermionic degrees of freedom and is diagonalized by 

fcir„,. „ _!_,■«. ^tl ,„;fU ^„„ffl^:„„fo „,. - Afc-(l+Acosfc) 



means of the transformation r/^ = -4= ^^ e''^' [afcC;-|-i/3fcc[] with coefficients 



a/c 



e Hamiltonian then has the form 
H = Y,kk^l^k-\Y.t^k■, Afe = Y/(l + Acosfc)VAVsin2fc. (2) 



and /3fc = —. — , ^Asmfc ^ r^^^ Hamiltonian then has the form 

^2[A|-(l+Acosfc)Afc] 



A. Correlation functions 

As specified within the next section, one- and two site- entanglement measures are obtained 
from the (one- and two-body) reduced density matrix whose entries can be related to various spin 
correlation functions 

M^w = {^stim) ■ g^^it) = {^\snt)siitm . (s) 

These can be recast in the form of PfaffiansSSiSLSS. Correlators defined in Eq. ^ have been calcu- 
lated for this class of models in the case of thermal equilibrium^^iSSiSSiSJ . In this case, the expression 



for the correlators reduces to the calculation of Toeplitz determinants. This is not the case here, 
since the initial state jf/;) is not an eigenstate of the Hamiltonian Eq.QJ. The correlation functions 
needed here have been obtained in22i2£. 
As initial state we consider two sites in a Bell state and all other spins being in the state | |) (or 

IT)) 

\n,)-^{4+e'''c^^\i^) (4) 

where the vacuum state is defined as |-IJ.) — |i . . . i). This initial state explicitly breaks the trans- 
lational invariance. The correlators i'^fASj^S^l'if,) can be expressed as a sum of Pfaffians as 
descibed inS*^ and allow, together with the magnetization, to evaluate the two-body reduced 
density matrix. Inaccessible by the above technique is the correlation {zL\S^\±). However, if 
(±1 Si'^{t = 0) |±) = 0, then it will remain zero during the subsequent evolution due to the parity 
symmetry of the Hamiltonian. We will consider exclusively this case in the present work. In the case 
7 = the particle number conservation leads to a considerable simplification since the vacuum is 
an eigenstate and the correlation functions can be evaluated directly without resorting on Pfaffians. 
For an open-ended isotropic chain dynamic correlators have been studied employing the pfaffian 
expression^. 



III. MEASURES OF ENTANGLEMENT 

On a qualitative basis entanglement is well understood. Both for distinguishable particles (e.g. 
spins on a lattice) and identical particles (e.g. free fermions/bosons).22. If a many-spin system is in 
a pure state, a rough measure for the entanglement between two subsystems is the mixedness of the 
reduced density matrix of the subsystem. We analyze the case when this subsystem is a single site 
and choose the one-tangle, which on site j is given by 

T(i)[p(i)]:=4detp(i) = i-(5;)' . 

Bipartite entanglement is encoded in the two-qubit reduced density matrix p^^^ , obtained from the 
wave- function of the state after all the spins except those at positions i and j have been traced out. 
The resulting p^'^' represents a mixed state of a bipartite system for which a good deal of work has 
been devoted to quantify its entanglemen1iSi24i2^. As a measure for mixed states of two qubits, we 
use the concurrence^ C[p'^'] 

C[p(^^] := max{0, 2Xjnax - trVR} ; R := p^^^Oy (g) Oyp^^'^ *Gy ® Oy (5) 

where \max is the largest eigenvalue of the matrix \jR,. For pure two-qubit states we have t*-^' = 
C^ =: t'^'. The concurrence, expressed in terms of spin correlation functions, is^^ 




a™ = 2max<; V(.g- tCJ' + (C ± O' - \/(J ^.gf^J^ - i(Mf tM^)2,o| (6) 

In the isotropic case, additional constants of the motion simplify the expression for the concur- 
renceSS. One- and two- site entanglement do not furnish a complete characterization of the en- 
tanglement present in spin chains. Following a conjecture put forward by Coffman, Kundu, and 
Wootters as follows 

5]cV<4detpW=rW. (7) 

it is interesting to study the difference of the quantities in Eq. (Q, interpreted as "residual tangle", 
i.e. entanglement not stored in two-qubit entanglement. This inequality was proved in Ref.^^ for 
a three qubit system, giving rise to the definition of the three-tangle as a measure of three-qubit 
entanglement. 
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FIG. 1: Top panel: The nearest neighbor concurrence C\ for the values of 7 = 0.1, 0.5, and 1 (from left 
to right). At A = 0.5 and A = 1, Ci has roughly the same value: 0.0264 and 0.0337 for 7 = 0.1, 0.1204 
and 0.1285 for 7 = 0.5; 0.2074 and 0.1946 for 7=1, respectively. For the Ising model at A = 0.9 we have 
Ci = 0.24. These values remain as uniform background in the concurrence of the singlet-type perturbation 
of the ground state (see bottom panel of Fig. O . Bottom, panel: The behavior of the separable point for 
generic 7 for left: 7 — > and right: 7^1. 



IV. GROUND STATE ENTANGLEMENT 



This section focuses on the ground state entanglement. In Fig. Q] upper panel, the nearest 
neighbor concurrence for the ground state is shown for different values of 7 as a function of the 
reduced coupling A. All exhibit a logarithmic divergence of the first derivative respect to A at the 
quantum critical point Ac = 1, and fit within finite size scaling theory. ^"^ This observation also 
applies to the first non-zero derivative of the concurrence at larger distance. The cusp, where Ci 
vanishes is the point where the large eigenvalues of i?, Eq. Q , for the invariant sectors (due to the 
parity symmetry) \Mz\ = and \Mz\ = 1/2 are equal. This indicates a change in the type of Bell 
state responsible for the entanglement. The value Aq where this degeneracy occurs (and the 2-qubit 
reduced density matrix is separable) converges from above to the critical coupling Ac = 1 for 7 — > 
as Ao = 1 -I- 7^-^^ and diverges as 7 ^ I as Aq = (1 — 7)"°-'^^ (see Fig. ^|. In the ordered region, 
the thermic ground state is a mixed state and as a consequence, the here shown concurrence is only 
an upper bound for the thermic ground state with spontaneously broken symmetry. However, it 
has been shown recentlji^ that if it is the triplet sector that furnishes the largest eigenvalue of i?, 
the concurrence is not affected by the thermic mixing. Therefore, the symmetry breaking does not 
affect the concurrence for I < A < Aq. 

The one-tangle for the ground state is shown on the left of figure |21 At the critical point for 
7 G [0.1 , 1] we could demonstrate that the sum of the two-tangles is much below the one-tangle. 
This is seen in the right plot of figure [21 If the CKW conjecture holds, this indicated that far 
the most entanglement is stored in higher tangles (yet to be quantified). This indication recently 
found further support coming from other indicators for higher entanglement as the scaling of the 
von Neumann entropy of a compact block of spins with the size of the block^^ and the amount of 
localizable entanglement, quantified by the maximum two-point correlation function^. 
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FIG. 2: The global tangle for the groundstate (left). For 0.1 < 7 < 1 we could verify the CKW conjecture 
for critical coupling (right). The one-tangle (thick line) is much larger than X^J^_^ C„ = 2 X^JJ^^ C„ (thick 
dash-dotted line) suggesting that the major part of entanglement should be stored in higher than two-qubit 
entanglement. The thin line is a guide to the eye showing that the sum of the two-tangles linearly tends to 
zero as 7 — > 0. 



DYNAMICAL EVOLUTION OF AN E-BIT 



The key to the time evolution of the original spin operators is that of the spinless ferniions c (t) = 
'^i[ai-j{t)ci — bi-j{t)cl], where the new coefficients are ax{t) = -^X^fe^osfca; (e*'^''* — 2z/3^sinAfct) 
and 

bx{t) = -j=J2k^'^^^^ '^kPkSinAkt. In the limit 7 = the previous expressions simplify consid- 
erably. In this case the magnetization, i.e. the z-component of the total spin S^ = '^^Si^, is a 
conserved quantity. In terms of ferniions this corresponds to the conservation of the total num- 
ber of particles, N = J2j''^3 = J2j'4^3' "'^°'' 7 ~* and |A| < 1 we find that ak — > and 
Pk -^ signfc. The energy spectrum is Ak = \l + \cosk\ and the eigenstates are plane waves 

1 ■ 



c, it) = 7^ E. E/ cos k{l - j>-'A^*q and 4 ^ -^Ei ^-"'c, 



Vn ' 



A. 7 = 0: The isotropic model 

In this section, we describe the dynamics of entanglement for 7 = 0. The Hamiltonian of Eq.lJQ) 
is then reduced to the XY model. Only in this case the ^-component of the total spin, S^ , is 
conserved. Consequently the Jordan- Wigner transformed fermionic Hamiltonian becomes a tight 
binding model for each sector with fixed S^. 



Propagation of states in the singlet sector 



We first consider the case of a chain initially prepared in a maximally entangled singlet-like 



state l^f,) on sites i and j as defined in Eq. 



Q. The state vector at later times is |^f,) 



E/W;(t)C; l-H) with Wi(t) 



/2N 



Ekie 



'-(t-i) 



cl-Vf 



Hj-1) 



jgiAfct 29 Pqj, g^j-^ infinite system, the 



coefficients become (up to a global phase) wi{t) — A={Ji^i{Xt) + e^'^ (~*) Jj~ii^t)}j where J„(a;) 
is the Bessel function of order n. For this initial singlet-like state the concurrence of sites n and m 
is Cn,m{i) = 2|w„(t)w*j(i)|. The time scale is set up by the interaction strength: the information 
exchange or "entanglement propagation" takes the time t ~ d/\ for d lattice spacings; i.e. the 
speed of propagation is A. The propagation of an EPR-like pair is demonstrated in Fig|21 where the 
concurrence is shown between two sites symmetrically displaced with respect to the initial excitation 
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FIG. 3: Left: Concurrence between sites n — —r, m — r, symmetrically displaced from their initial position 
i = — 1 and j = 1 {ip = n). Right: Summed concurrences for an initially entangled site (n — 0, full line) and 
an initially unentangled site (n = 10, dashed line). The initial state is a singlet {(p — n) created on sites 
1 = and j = 1. 

(Cr,-r(i) for i = —j = 1) together with the sum of concurrence Ctot.n = X^m ^n,m of the single site 
n. The latter tends to the same stationary value (around 1.6) for arbitrary n, independently of the 
initial conditionsSS. This means that the initial state becomes homogeneously spread at long times 
and so does the concurrence. 

The same structure as for the concurrence is also found for the one-tangle. This fits within the 
picture suggested by the CKW conjecture^'^. For the isotropic model it is possible to analytically 



check this conjecture. We obtain 4detp^- ' = 4|wjp(l - \wj\'^) and X^n^^i '^],n = 4l]„ 



4|w;ji2 



^/i 



\w-,w„ 



(1- 



). This has already been observed for these W-type statesS^. Thus, the two quantities 



coincide, meaning that the entanglement present in the system is restricted to the class of pairwise 
entanglement, and no higher order entanglement is created. Therefore, the information in the 
one-tangle is already contained in the concurrence. 

In order to gain information on what type of Bell state propagates and on eventual state mutation, 



we study how similar is the mixed state p, 
the fidelity 



(2) 



to the initial \^!^). This similarity is quantified by 



Fn.At) = Tr{p(2)^(t) 1*'^) (vl/^l} = 1 h«„(t) + e 



-l(p 



Wj 



M 



(8) 



It was shown to display in-phase oscillation with respect to one-tangle and concurrence^^. Thus, 

(2) 

when the entanglement wave arrives, pn,m becomes more similar to the initially prepared state. 
That means that the state itself is propagating along the chain, taking with it its entanglement. 
This propagation is far from being a perfect transmission, due to the entanglement sharing with 
many sites at a time. 



2. Propagation of states in the triplet sector 



We next analyzed the propagation of the state 






1 



^^7=2^' 



e'''clc])m 



(9) 



These are not single-particle states and since they are superpositions of components pertaining to 
different spin sectors, one cannot take full advantage of the conservation of the magnetization. As a 
result, the concurrence is given by C = max{0, C^^^ , C^^^}; a concurrence of the form C^^^ indicates 
that |\I/)-like correlation arises, while for C^^' > C^^-' the correlations between the two selected sites 
is more |$)-like. This already suggests that changes in the propagated type of entanglement could 
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FIG. 4: Upper panel: Time evolution of the concurrence between sites —n and n for the initial state 
I $'^5 g). The plot is cut at 0.2 in order to make the revival after the crossing visible. Lower panel: Fidelity 

(whiteihigh, blackdow) of the states |<I)'^°p*) (left) and |\l>'^°p') (right) in the density matrix p_^ ^(t) for the 
same initial condition as above. 



occur. 

The entanglement still propagates with the sound velocity ~ A along the chain, as can be seen 
in Fig. Q, where C-n,n is displayed for an initial state |$f.5 5). At the intersection between the 
sites of the initial Bell state a revival takes place, and the entanglement then also spreads out from 
this intersection point. In contrast to the previous case, when the two initially entangled sites are 
separated by an odd number of spins, the propagating quantum correlations change character. This 
is seen from the fidelity of the Bell states \^'^ „) and |^^j „) in the time evolution of the initial Bell 
state l^f ,). These two fidelities are maximized by a proper choice of the phases^S. which depend 
only on the initial and final positions i, j, m and n and are shown in the bottom of Fig. Q ~ left 



plot: l^vopty^ right plot: l*"^" 



for i = —j = 5 and n = —m. Since the fidelity of |$* 



the 



vacuum is 0.5 due to its component |ii), it is seen that after the crossing point the state has become 
|$)-like. In fact, at the single-site crossing, the amplitude for two parallel spins cannot survive and 
the outgoing states of this scattering event only contain antiparallel spins. If the initially entangled 
spins are separated by an even number of sites, the crossing involves two sites and the character of 
the state is preserved. 



B. 7 7^ 0: Singlet onto the vacuum 



We now study the model for generic 7 and infinite chain. In contrast to the previous section, here 
the complexity of the computation grows with the distance d of the sites because the dimension 
of the Pfaffian expression for the correlation functions is 2d x 2d. For the critical Ising model it 
turns out to be sufficient considering d = 1, since the concurrence vanishes for larger distances. The 



initial state is the singlet 1*1^2) j Eq. 0). 



Concurrence 



Figure shows the nearest neighbor concurrence for 7 — 0.5, 1.0 (1st and 2nd panel) and 
A = 0.5, 1.0 (left and right, respectively). A rough estimate of the propagation velocity can be 
taken from the contour lines in the plot. It coincides with the sound velocity, which is roughly A as 
for the isotropic model. For increasing 7 it slightly decreases for eventually returning to ^s^SSiSiiS. 
The concurrence on the original singlet position decays quickly and the damping of oscillations is 
stronger for increasing A. We notice an instantaneous signal which, sufficiently far away from the 
initial singlet position is spatially uniform. This phenomenon reflects the creation of entanglement 
from the vacuum, which is characteristic for the anisotropic XY models. It is originated from 
the double spin flip operators 7A/2 ^^ ^t^''i+i (a-nd its Hermitean conjugate), which are absent in 
the isotropic model. The initial slope and the type of Bell state created agree well with simple 
perturbative arguments. This pure vacuum signal dies out very quickly. Towards the critical 
coupling and the Ising model, the damping of the concurrence propagation gets stronger. The 
vacuum signal survives much longer for medium A and 7^1 such that for 7 = 0.5 and 7 = 1 it 
interferes with the propagating singlet. Although the damping of the propagation becomes stronger 
at the critical coupling, nevertheless it is of pure dynamic origin and not related to the quantum 
phase transition, since the initial state is not the ground state (where the critical behaviour is 
encoded in). Consistently, the damping turned out to be independent of the size of the chain. 
The concurrence for 7 = 0.5 is shown in the top pannel of Fig. [5] we see a clear propagation of 
the concurrence, which is only slightly stronger damped than for the isotropic model, but there is 
a creation of entanglement from the vacuum. A "shoulder" appears in the singlet peak of Ci^i+i, 
i.e. on the original singlet position, which is due to triplet-type entanglement^^. The squares of the 
preconcurrences C*-^' and C'-^\ shown in the 3rd panel of Fig.[Sl alike the analysis of the fidelity of 
the Bell states in the propagating signalSSi, demonstrate that the propagation is of the same type as 
the initially created Bell state (here in the singlet sector) , whereas the concurrence created from the 
vacuum is of triplet type. For the transverse Ising model (2nd panel in Fig. [SJ, the shoulder on the 
original singlet position and the vacuum creation gets even more pronounced, but all the signals die 
out much quicker. For A = 0.5 one cannot speak any more of a clearly propagating entanglement 
signal. At the critical coupling, the propagating signal almost completely disappeared. The next- 
nearest neighbor and next-next-nearest neighbor concurrence are shown in Fig. |B1 Both show a 
narrow wall created from the vacuum, which for C3 :— Cx-3,x is broader. C3 unveils an additional 
feature: it shows a large contribution at t around 4 at a; = 3. It indicates an EPR-type propagation, 
which we observed already for the isotropic model (see Fig. O. Considering the velocity A, the time 
of appearance indicates that the two fragments have to cross first. At and near the critical coupling 
C2 and C3 identically vanish on the domain of the demonstrated plots. 



2. The global tangle 

Is it possible to understand how the entanglement, originally stored into the singlet, is going to 
share among all the spins in the chain? In order to get an idea about what might happen, one can 
study the one-tangle. For small anisotropy parameter, this quantity is qualitatively very similar 
to the concurrence: there is an entanglement wave, propagating with velocity A. For growing 7, 
the wave is suppressed in favor of a homogeneous growth of the one-tangle which saturates on very 
short timescales. This signal cannot be due to the initial singlet, of course, but it is created from the 
vacuum. This was confirmed by looking at the total entanglement created solely from the vacuum, 
i.e. the initial state being the vacuums^ The right-most plot in Fig[7| shows the residual tangle for 
the Ising model at A = 1; there is no violation of the CKW conjecture. 
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FIG. 5: Upper two panels: The nearest neighbor concurrence for medium A = 0.5 (left) anisotropy 7 — 0.5 
and A — 0.5 and 1 (left and right, respectively). The propagation is stronger damped than for small or 
zero anisotropy. There is a creation from the vacuum which reaches the propagating pulse. At the critical 
coupling (right), the nearest neighbor concurrence from the initial singlet dies out immediately and so does 
the vacuum creation. For the Ising model, only few bumps are residues of a propagation . 3'rd panel: The 
square of C'^' (left) and C'^' (right) for the Ising model at A = 0.5 demonstrate that the propagating 
signal is in the singlet sector, whereas the background is in the triplet sector. Bottom panel: Cx-i,x for 
the perturbed ground state and the Ising model: A = 0.5 (left) and A = 0.9 (right). There is a background 
concurrence corresponding to the ground state value (see Fig. 0. A propagating signal is seen in a valley 
of extinction. 
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FIG. 6: Left - The next nearest neighbor concurrence C2 for the Ising model far from the critical coupling. 
It is zero at and close to the critical coupling. Right - The next-next nearest neighbor concurrence C3 for 
the Ising model far from the critical coupling. Also C3 vanishes at and close to the critical coupling. At 
A = 0.5 we see a considerably large signal, which is due to an EPR-type propagation of a "split" singlet. 
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FIG. 7: The total tangle of site x for "/ — 1. In contrast to the nearest neighbor concurrence there is a clear 
propagating signal here fro A = 0.5 (left). It is mounted on top of a non-zero background signal coming 
from the vacuum. At critical coupling (middle) the vacuum background of the one-tangle grows up to about 
0.75; the propagation is hardly visible. Right: The CKW conjecture applies: the "residual tangle" indeed 
is positive. 



Next, we will study the relative deviations in the one-tangle respect to the time-evolved vacuum 



At: 



rel 



det 



P( 



(1) 



detp 



(1) 



(10) 



We want to stress that the time evolution operator does not in general preserve the relative order 
induced by an entanglement measure on the Hilbcrt space. This manifests in negative values of 
Atj . Eventually, this is due to the fact that superposing (as well as mixing) orthogonal maximally 
entangled states of the same type diminishes the entanglement. This may lead to a negative Ar'''^'. 
It is worth noticing that the singlet state is inserted into the vacuum; this is not a superposition 
of the vacuum and some other state. We choose the same parameter range as above. Along the 
axis t = and A = and for site numbers larger than two (and smaller than one) we have that 



.(1) 



,(1) _ 



detp^ "^ = 0. In these cases, also det/9L^^^ — 0, and we chose the plotted value being zero in these 

cases. The analysis of Ar'"'^' tells us that for small anisotropy and sufficiently far from the critical 
coupling the global tangle is dominated by the local perturbation of the vacuum by the singlet 
(Fig.lHl top), meaning that the total tangle is concurrence dominated. For the isotropic XY model, 
the global tangle was given entirely by the sum of the 2-tangles such that the CKW-conjccture 
would conclude that there is no higher tangle contained in the system. In the presence of a small 
anisotropy, this is no longer true, in particular near to the critical coupling A = 1 (Fig. (SJ. Rising 
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FIG. 8: Upper panel: The relative tangle deviation Arg "^ from the Vacuum tangle at site number 6 for 
different values of A for a fixed value of 7: 7 = O.f (left), 7 = 0.5, and 7 = 1 (right). It is plotted as a 
function of the reduced time r = t/X; it is nicely seen that the oscillation frequency grows linearly with 
A. Lower panel: The same as above but for the ground state and 7 = 0.1 (left) and 7=1 (center). The 
oscillation frequency grows linearly with A. Right: anomaly at t=0: it demonstrates the non-local impact 
of the singlet-type perturbation at the critical coupling. 



anisotropy and A — > 1 enhance the vacuum domination. 



C. 7 7^ 0: Singlet-type perturbation of the groundstate 

One could wonder, whether or not the results in the preceding sections were specific to the vacuum 
state or what would happen for different states. We therefore discuss in this section the propagation 
of a singlet-like perturbation of the groundstate \GS), i.e. the time evolution of the initial state 



\S)^:^—ic^-C2)\GS) 



(11) 



We note that this state differs from a singlet on the ground state since the operators Ci create global 
excitations. It turned out that the resulting state is nonetheless very similar to the ground state 
itself. For the isotropic model \GS) = Ifl") and hence the dynamics is the same as on the vacuum. 
Apart from the propagation of the pulse with velocity about A, there are several qualitative 
differences for the concurrence. One of them is that the propagating pulse and the initial Bell state 
on sites 1 and 2 is the triplet with zero magnetization. This comes from subtleties in the Jordan- 
Wigner transformation^^. The concurrence plots in the bottom panel of Fig. |Slshow that the inital 
state is almost fully entangled on sites 1 and 2. This demonstrates how close is the ground state to 
|-f|-). Two new features appear in the concurrence signal. Firstly, the entanglement pulse propagates 
on top of a nonzero background level, which in very good agreement coincides with the nearest 
neighbor concurrence of the ground state, being around 0.2 at both, 7 = 0.9 and 7 = 0.5i^ (see 
Fig. nj. This shows that not only is the ground state very close to the vacuum, but that also the 
state \S) is very similar to the ground state as far as nearest neighbor concurrence is considered. 
Secondly, in contrast to the singlet on the vacuum, we here have to deal with a propagation that 
eliminates the background concurrence. The latter feature gets more pronounced when approaching 
the Ising model and critical coupling (see bottom panel of Fig.O. This elimination is understood 
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from the fact that joining entanglement of the same type (here: two-site entanglement) in form 
of states that are orthogonal to those forming the entanglement already present, diminishes the 
entanglement. Whereas here the impact of the propagating concurrence pulse coming from the 
initial 0-triplet is the stronger the closer we get to the critical coupling and the quantum Ising 
model (the opposite of what we observed for the initial singlet onto the vacuum), the situation 
would be the same as far as the total signal along the propagation line is concerned. This in fact 
gets more suppressed with growing A and 7. 

For the global tangle we find essentially the same behavior as for the singlet on the vacuum. Only at 
very short times the singlet-type perturbation deminishes the global tangle. In the bottom panel of 
FiglSlwe choose two different anisotropics 7 = 0.1 and 1 and compare AtJ'^'' for different couplings A. 
The short-time behavior of Arg "^^ shows a marked anomaly at the critical coupling. This is a cursor 
of the non-local impact of the singlet-type perturbation of the ground state, which is pronounced 
at Ac = 1. At 7 = all curves should eventually tend to zero. 



VI. CONCLUSIONS 

In the present work we studied the effect of a singlet-type perturbation on the entanglement of 
an infinite spin system. We considered quantum XY models for general anisotropy. The dynamics 
of entanglement was studied as function of the distance to the local perturbation at i = 0, and 
of the reduced interaction strength A (up to the common quantum critical point of the models at 
A — 1). For this class of models we analized a conjecture formulated by Coffman, Kundu, and 
Wootters quantifying the weight of the pairwise relatively to the global entanglement, measured by 
4 det pi . In all cases the main propagating signal is in the same sector as that one initially created 
and the propagation velocity is in good agreement with the sound velocity of the model, which is 
roughly A. The isotropic model, i.e. anisotropy 7 = 0, can be mapped onto a tight binding model, 
and consequently entanglement propagates only, remaining pairwise. In addition we have an EPR- 
type propagation of the concurrence. We found that states in the triplet sector (with concurrence 
C = C*-^-*) are mutated after a crossing on a single site. The global tangle and the concurrence (whose 
square is the 2-tangle) satisfy the Coffman-Kundu- Wootters conjecture with zero residual tangle. 
This means that the system contained only pairwise entanglement, measured by the concurrence. 
The Hamiltonian does not create any entanglement; it distributes the initially created pairwise 
entanglement. Also for general anisotropy we found evidence for an EPR-type propagation. The 
propagation is suppressed, compared to the isotropic model. The suppression is the stronger the 
closer the system is to the critical coupling and the quantum Ising model. For the latter we found 
a very rapid damping of the singlet in the nearest neighbor concurrence. For all larger distances 
we considered (up to 7 lattice spacings) the concurrence is zero. A peculiarity of the anisotropic 
models is the instantaneous creation of concurrence from the vacuum all over the chain. It decays 
very quickly when approaching the quantum critical coupling and getting close to the Ising model. 
Neither effect is of critical origin, though. Comparing quantitatively 4 det pi and concurrence we 
conclude that for the anisotropic model the propagation of the concurrence is a small effect respect 
to the creation of higher tangles. Medium interaction strengths and/or small anisotropy favor the 
propagation of the singlet. For the singlet-type perturbation on the ground state, the concurrence 
signal occurs along a "valley" in the constant background concurrence. The propagating extinction 
gets more enhanced with growing anisotropy and approaching the critical coupling; for the critical 
Ising model a weak propagating signal in the valey remains. The background nearest neighbor 
concurrence coincides with that of the ground state, indicating that the nearest neighbor concurrence 
away from the initial perturbation is unaffected, and that the Hamiltonian cannot notably create 
nearest neighbor concurrence beyond that level. The dynamics of the total tangle is basically 
unchanged. In the short-time behavior of Ar^^' an anomaly at A = 1 unveils the non-locality of the 
singlet-type perturbation of the ground state. From the perspective of quantum information, our 
results can be read as a transfer of a unit of pairwise entanglement (an e-bit) within a spin chain 
by Hamiltonian action. The so transported concurrence can then be distilled at the destination 
poinliiS^. 
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